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The more exact upper estimate of the percolation threshold for the 
site problem on the quadratic lattice Z 2 have been found on the basis of 
'. the cluster decomposition. It is done by the number estimate of cycles 

on Z 2 which maybe external boundaries of finite clusters. 
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Generally, percolation theory studies the connectedness relation for 
random sets in topological spaces. Problems representing greatest in- 
terest in the theory appear in noncompact topological spaces. 

The existence of a noncompact connected component in random 
set realizations with nonzero probability is the crucial question in the 
theory [1]. In such a general setting of the problem, however, it is im- 
possible to obtain any strong results. That is why the basic objects of 
the percolation theory are random sets in M. d and Z d at d = 2, 3 gen- 
erated by stationary random fields [1, 2]. The first case corresponds to 
so-called continuous percolation theory and the second one - to discrete 
Q_c percolation theory. This limitation is a consequence of the fact that it 

^ r-| is the problem set in such a way is mainly needed in physical applica- 

tions. Therefore, it may be considered as the object of mathematical 
physics. However, despite the above restriction, the basic problem is 
very complicated and hardly lends itself to exact mathematical study. 
Rigorous results referring to discrete percolation theory which is its 
mostly developed direction have been summarized in the monograph 
[3]. Later results and also results referring to continuous theory have 
been summarized in the excellent review [4] . Physical literature dealing 
with different nonrigorous heuristic approaches to the investigation of 
the percolation phenomenon, computer experiments and special appli- 
cations of the theory in theoretical physics is giant and we shall not 
pay attention to it in this short introduction to the problem. 

Despite of the great flow of publications at 80-th in physical litera- 
ture and simultaneously appearance of some reviews in which mathe- 
matical results are summarized, main questions of the percolation the- 
ory considered as the object of mathematical physics remain without 
answers up to now. It is the situation even for simplest case being 
studied in the theory when random set is generated by Bernuolli's field 
on Z 2 . Even for this case any algorithm of the percolation threshold 
calculation (see §2) and, more generally, any algorithm of percolation 
probability evaluation was not found for an arbitrary periodic graph [1]. 
Our communication is dealing with the revision of upper estimate of 
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percolation threshold as compared to what may be found in literature 
[4, 5] for the periodic graph on Z 2 which is called the square lattice. 

2. Percolation theory problem on Z 2 

Let us consider an infinite graph with the vertex set Z 2 . For sim- 
plicity of statements and arguments, we shall study it as one immersed 
into R 2 . The adjacency relation on the graph is defined by the set 
of pairs {(x, y) G Z 2 : xipy} where Xipy if and only if y = x ± e! or 
y = x ± e 2 , ei = (1, 0), e 2 = (0, 1). We shall follow the terminology of 
statistical physics and that is why we shall call such a graph the square 
lattice and denote it by the same symbol Z 2 . 

Let {c(x)} is the Bernoulli random field with the parameter 

c = Pr{5(x) = 1} 

which is said to be the concentration. Below, the tilde marks the ran- 
domness of used objects. The field {c(x)} induces the random set with 
realizations {M : M C Z 2 } where M = {x : c(x) = 1}. Sometimes 
we shall call these realizations configurations of completed vertices or, 
simply, configurations. This probability distribution for possible real- 
izations M is clearly determined completely by the following probability 
collection 

Pr{M : A C M} = c |A| , AcZ 2 . 

Here and further | • | = Card{-}. 

Naturally, the adjacency relation if induces the connectedness rela- 
tion for vertices having included into any configuration M. We shall 
call two vertices x and y linked in M if there exists a path (xf, % = 
0, 1, 2, n), Xi G M, x = x,x n = y and Xiipx i+ i, i — 0, 1, 2, n — 1. 
The connectedness of vertices of M is the equivalence relation. There- 
fore, each configuration M is uniquely decomposed on some disjoint 
equivalence classes M — generated by the connectedness rela- 

tion of vertices. These classes we shall name clusters. The collection 
of clusters related to configurations M will be denoted by W[M] = 
{Wj]j G N}. If x G Wj for a j G N in configuration M, then we shall 
denote this cluster Wj by W(x) [3]. 

Let us introduce the random field {a(x); x G Z 2 } on the base of the 
set M, 




1; x G W C M, \W(x)\ < oo, 
; otherwise . 



In view that {c(x)} is the uniform field (i.e. the probability measure is 
invariant relative to translations on vectors (niei + n 2 e 2 ), ni,n 2 G Z), 
the random field {a(x)} is uniform too. Therefore, the probability 

Q{c) = Pi{a(x) = 1} 
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does not depend onx e Z 2 . This probability is a nondecreasing function 
of c [3]. If Q(c) > 0, then the percolation on the {c(rc)} is said to 
exist. In connection with this fact, the following characteristic value is 
introduced 

c* = sup{c : Q(c) = 0} . 
c* is called the percolation threshold. 

3. Finite clusters on Z 2 

At construction of some upper estimates of the probability Q(c), it 
is necessary to enumerate all finite clusters W{x) including the fixed 
vertex x. We perform the enumeration by using the concept of the 
external boundary for each finite cluster. To this end, let us introduce 
on Z 2 the following new concept of the adjacency relation <p side by side 

2 

with the adjacency relation ip [3,4]. We denote Z the corresponding 
graph having the same vertex set but with the adjacency relation (p. 

Vertices x and y are named the ^-adjacent ones if one of two cases 
takes place: 1) xipy, 2) either y = x + ei ± e 2 or y = x — ei ± e 2 . 

The relation <p induces the new relation of vertex connectedness on 
each configuration M. This connectedness relation is the equivalence 
relation too and it leads to a decomposition of M on some connected 
sets of vertices relative to (p. 

Let us introduce the external boundary concept of a finite cluster 
on Z 2 . 

Definitionl. The set dW is named the boundary of the cluster 
W if it consists of those vertices y which are <p-adjacent to vertex x in 
the cluster W but do not belong to it. 

Definition2. The external boundary dW of the cluster W is 
the set of vertices u G dW such that for each of them there exists an 
infinite (p-path a(u) on completed vertices ofl? and, moreover, u is the 
unique vertex in a(u) belonging to the union W\JdW . 

Further classification of finite clusters W(x) is performed by the enu- 
meration of all possible external boundaries dW(x). Following state- 
ment is the key for such enumeration. It repeats the corresponding 
statement in the monograph [3] with the exception of the last item. 

Theorem 1. Let W{x) be a finite cluster for a fixed vertex x, 
\W{x)\ < oo. Then W(x) has a nonempty finite external boundary 
dW(x) having following properties. 

1. dW(x) is the (p-connected vertex set in Z 2 which presents the 
cycle i.e. dW(x) = (xi, X2, x n ) where X{ ^ Xj, i ^ j, n = 
\dW(x)\, Xi<pXi+i, i = 1,2, ...,n, x n+ \ = xq and each vertex has only 
two ip-adjacent vertices in the 7. It is possible to introduce the definite 
orientation on dW(x). (it will be used the counter-clockwise orientation 
in further arguments.) 



3 



2. The vertex x is contained in the finite set lnt[dW(x)] defined by 
lnt[dW(x)} = {u <£ dW{x) : V(a(«) : \a(u)\ = oo) (a(u) n dW{x) ^ 0)} ; 

3. Let u,v,w G be three (p-adjacent vertices following one 
after another according to the introduced orientation. Then ifwpv, the 
vertex w belongs necessarily to one of the following collection: 

a) in the case uipv the set of three elements {2v — u, 2v — u ± e} 
where e is a unit basis vector that is orthogonal to the vector (v — u); 

b) in the case v — u + e(ei — e 2 ) the set of five elements 

{2v — u, 2v — u + eei, 2v — u + ee 2 , 2v — u ± e'/2} 

m i/ie case v = u + e(ei + e 2 ) ; 

{2t> — u, 2v — u + eei, 2f — w — ee 2 , 2v — u ± e'/2} 

where e' &e an orthogonal vector to (v — u) and |e'| = \/2 and e = ±1. 

The proof of first and second statements of Theorem 1 is obvious (see 
[3]). What about the last statement, it is easy to verify its justification 
from the explaining figure l(a, b). The full rigorous proof is found by 
us; it is quite tedious and goes beyond the scope of this communication. 
Unlike the proof of the analogous theorem in [3] , our arguments are not 
used an application of Jordan's theorem. Further we omits the proof 
and pass to the proof of the main statement in next paragraphs. 

4. The cluster decomposition on Z 2 

We shall consider the probability 1 — Q(c) = Pr{a(0) = 0}. Let 
A = {W : W = W(0) is cluster, \W\ < 00} be the collection of all 
finite clusters including the vertex 0. We define the event 

A{W) = {M : W = W, G M, W e {WfJ G N}} 

for any cluster W G A. This event has the definite probability 

Pr{A{W)} = c\ w \l-c)\ dw \ . (1) 

According to statements in the previous section, any cluster of this 
collection A is corresponded to a cycle 7. Vertices of this cycle are (p- 
adjacent and such that G Intfy]. In this connection, let us introduce 
in consideration the collection B of all y>cycles having the last property. 
Let 5(7) be the event 

5(7) = {M : G M,W(0) G W[M],dW(0) =7} (2) 
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defined for any </>cycle 7 G B. It is represented by the finite union of 
mutually disjoint events 



B{l) = |J A{W). (3) 

weA:dw=-r 

Hence, according to (1) and (3), such event has the definite probability 

P( 7 ) = Pr{5( 7 )}. 

It is equal 

P( 7 ) = ^ Pr{A(Vy)} = ^ c |H/| (l-c)W. 
Let us note that 

{5(0) = 0} = |J A(W) . (4) 

W&A 

The collection A is decomposed on some disjoint classes of clusters. 
Each class consists of those clusters W E A which have the same exter- 
nal boundary. It is realized dW = 7 for them. Therefore, the following 
representation is correct 



U ... = u U ■■■■ 

WeA iet3 WeA:dW=-y 

Then we obtain 

{ a(x) = 0} = |J Bin) 

using (3). Finally, we come to the statement 

Theorem 2. The probability 1 — Q(c) is represented by the decom- 
position 

l-Q(c) = £P( 7 ). (5) 

7GB 

Usually, such decomposition is said to be the cluster one in the 
percolation theory [4] . 



5. The main theorem 



The cluster decomposition (5) is represented by the sum of prob- 
abilities of some disjoint events. Therefore, definitely, the cluster de- 
composition is convergent. The function Q(c) is not equal to zero only 
if c > c* > 0. So, it is not an analytic function on the concentration 
c and c = c* is its the singular point. Earlier [4], it was obtained that 
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6/7>*>l/3. We give the improvement of the upper estimate at the 
below-formulated statement. 

Theorem 3. The inequality c* < c = 3 — y/E is correct for 
Bernoulli's random field on I? . 

Proof. We use the elementary estimate 



P( 7 ) < (1 - c) 



hi 



that follows from the Definition 2 and Eq.(2). Using it and Eq.(5), we 
come to the upper boundary 

E p w < E^ 1 - c )' 7 ' = E ^ - c ) 2 ^ ( 6 ) 

7e6 yeB ~ t el3 -.\~i\=2k 

where = Card{7 G B : | 7 | = 2k}, k > 2. 

Further, we shall find the upper estimate for the value r^. Consider 
the infinite path a(0) = (je^ j = 0, 1, 2...) with initial vertex 0. Then, 
according to Theorem 1 (2), each cycle 7 G B necessarily crosses the 
path in some vertices. Let us select the vertex in this set of all inter- 
section vertices which is nearest to the vertex 0. Denote it by z T All 
cycles of B is decomposed on disjoint classes C 7 , i.e. the cycles having 
the same vertex z 7 belong to the same class. This class decomposition 
induces the decomposition of the cycle set {7 G B : | 7 | = 2k} on corre- 
sponding classes Cf^ where / is the distance from to z 1 . In addition, 
/ < k and therefore, 

{ 7 G£: | 7 | = 2A;} = [jC?\ r k = 

1=1 1=1 

Let 7 = (z y = x , xi, X2k-i, %2k — z -y)- According to the earlier 
introduced orientation, the vertex x\ that follows after x in the cycle 
7 may be one of the set only 

(z 1 + ei, Zj + ei + e 2 , z 1 + e 2 , z 1 - e x + e 2 ) . (7) 
Using this ordering, cycles belonging to the same class are dis- 

(k i) 

tributed on disjoint collections C\ , % = 1,2,3,4 independence of the 
selection of the vertex x 1 in the cycle 7. Hence, 



(fc) 

! 

1,2,3,4 i=l,2,3,4 



c<" = u <f>, cf = £ |c« 

i=l,2,3, 

We obtain as the result 



(k,i) 



Z=l 1=1,2,3,4 
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We must find upper estimate of the value 

Ak,i) 



C 



(k,i) 



It is easy to see that 



C v[ 2k where is the set of paths beginning at the 



vertex z 7 and the distance from z 1 to the vertex along the path a(0) 
is equal to /. Besides, they have the length 2k — 1 and the vertex x\ is 
the i-th in the finite sequence (7). We sort out paths 7 with the length 
2k — 1 but not with 2k. It is connected with the fact that the last edge 
X2k-i<fXo °f the cycle (xo, xi, X2k-i, #o) is fixed by the collection of 
vertices xq, x±, X2k-\- Therefore, 



C 



(k,i) 



< 



(2fc-l,i) 



S2fc-1 • 



It is easy to see that s n does not depend on / and i due to our con- 
struction. Let us represent 



s -- s r 



(8) 



where 



s+ = Card{7 G V\ n,i) : x n _^x n } , 
s* = Card{7 G V\ n,i) : :r„_iy? x n , x n ^(px n }. 

We introduce the two-component vector [s£,s£] (not distinguishing 
column- vectors and tuple- vectors) . Then the following equation takes 
place 



n-l 

x 

n-l 



with the transfer matrix 



1 2 

2 3 



according to Theorem 1 (3). Hence, 



-n-l 



(9) 



Eigenvalues of the matrix T is equal to A + = 2 + v^, A_ = 2 — \/E 
and corresponding eigenvectors (they are orthogonal but are not nor- 
malized) have the form 



t+ = 



1 

2 



1 

i-VE 



Decomposing the vector [si",s*] = [2,2] on eigenvectors t + ,t_, 

[*i">*i] = 9+t+ + g-t-, 
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we obtain g + — 1 + n/5/5, g_ = 1 — Therefore, [s£,s£] = 

(? + A" _1 t + + g_A_t_ and, according to (8), (9), we find 

s n = {g + \r 1 + g-\l- 1 ) + \[g + \r\l + ^)+g-\ n _-\l-^>)) < 

<4(2 + Vf) n l . (10) 

Since r k < 4(/c — l)s 2 fc-i, hence, putting n = 2k — 1 at the estimate 
having found, we obtain the majorant series from (6). It is summable 
at c > 3 — \/E. Then, if this condition is satisfied, we shall obtain 

£p( 7 ) < oo. (11) 

Completion of the proof is performed on the base of arguments which 
are standard in the percolation theory (see, for example, [5]). The in- 
equality (11) permits to apply the Borel-Cantelli lemma (see, for ex- 
ample, [6]) to the event family B. According to this statement, the 
probability of the event that consists in simultaneous realization of an 
infinite set of events belonging to the family, is equal to zero. Then, 
with the probability one, there exists a maximal cycle 7 G B. For 
this cycle, there is a vertex z out of it and, besides, there is the in- 
finite path a(z) without intersections and with the initial vertex z. 
The path does not intersect the cycle. This means that the event 

{M : 3 (w G W[M]) (|#| = 00)} of Bernoulli's field {c(x)} has the 

probability 1. On the other hand, the countable decomposition 

{M : 3 (We W[M]) (\W\ = 00)} = [J {M : W(v) G W[M], \W(v)\ = 00} 

(12) 

takes place. In view of the random field {a(x)} is uniform, the proba- 
bility 

Pr{M : W(v) G W[M], \W(v)\ =00} = Q(c) 

does not depend on v. Therefore, it cannot be equal to zero, since the 
following inequality is correct according to (12) 

1 < ^Pr{M : W(v) G W[M},\W(v)\ =00}. ■ 

Consequence. At > c , the probability Q(c) may be represented 
with any preassigned accuracy by finite sum of series (5). It is defined 
by summands which correspond to cycles 7 with the length not exceeding 
2m with fixed m G N. In this case, the following error estimate takes 
place 

E p w ^ E (i--) l7l <E( 1 - c ) 2fc " ls ^ = 

7GH:|7|>2m ~f£]3:\y\>2m k>m 
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oo 



-i 2fc-2 



16(1 -c) 2 J2 (k-l)[(2 + Vt){l-c) 



k=m+l 




< OO 



Thus, the decomposition (5) is the key for solving of the main perco- 
lation problem on Z 2 that is to calculate the probability Q(c) with any 
guaranteed accuracy. In addition, due to presence of the singular point 
= c*, the majorant series for the remainder must have a singularity in 
a point o > c* for any initial finite sum of the cluster decomposition. 
Such a situation will take place in any way of upper estimation of the 
percolation threshold. Then, it follows that the solution of the main 
problem is impossible without building of a calculation algorithm for 
the percolation threshold c* with any preassigned accuracy. 



One can see on the base of the above proof that the singularity * is 
possibly connected with the singularity of the generating function for 
the nonintersecting contours number. If it is true, so the regular method 
of sequential approximative calculation of this point may be based on 
the sequential exclusion of intersections in paths (zj,xi, ...,X2k-i) cor- 
responding to cycles 7. Such an exclusion must be fulfilled with the 
increasing of approximation order. In this case the estimate given by 
Theorem 2 may be considered as the zero approximation. 

Generally, the algorithm of approximations of the point c* should be 
consist in such a construction when estimates ~ < c* < c*, c~ +1 > c~, 
c+ +1 < c+ are build sequentially step by step. They have the property 
lim = c* and besides, for sure, Q(c) = at c < c~; Q(c) > at 



6. Discussion 



c> a 



,+ 

'n ■ 
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